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Abstract 



Via the numerical analysis on the intraband dynamics of optically excited 
semiconductor superlattices, we hnd that time-integrated squared THz emis- 
sion signals can be used for probing both dynamic localization and multi- 
photon resonance in the coherent regime. Competition effect between dy- 
namic localization and photon-assisted transport has also been discussed. 
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Recently, semiconductor superlattices (SL) driven by external electric fields have been a 
subject of intensive investigations both theoretically and experimentally. The most signif- 
icant example is Bloch oscillation (BO) which was theoretically predicted many years ago 
and experimentally demonstrated recently in SL driven by an uniform electric field. While, 
in the regime of tight-binding approximation, an intense laser field (ac fields) can make a 
charged particle localized around the initially injected place, which is the so-called dynamic 
localization Then the natural subsequent question will be the following: What will hap- 
pen if the electrons in SL is subject to the combined dc-ac fields? Actually, this question is 
answered by the recent several theoretical predictions and experimental observations. For 
example, multi-photon absorption 0, absolute negative conductivity |J, photon-assisted 
transport |3]], fractional Wannier-Stark ladders [g], etc. Particularly, the recent experiment 
on the multiple quantum-well superlattices driven by dc-ac fields has confirmed through 
probing the static I-V characteristics that the electrons in SL can tunnel through the adja- 
cent wells through the stimulated emission and absorption of photons [|[]. In this facet, the 
common action of the dc-ac field plays the constructive role on the transport properties in 
SL. However, the dc-ac field can also play the destructive role on the transport in SL, which 
can be manifested in dynamic localization induced by dc-ac fields, provided that the ratio 
of the Stark frequency uj b {= eF d, d is the SL lattice constant and F is the amplitude of 
dc fields) to the ac field frequency us is an integer n, and edFi/Tiu (Fx is the amplitude of 
ac fields) is a root of ordinary Bessel function of the order n: J n (edFi/huj) = 0. In this 
case, the electrons will become localized. Hence, the effect from both dc and ac fields play 
the dual role in the transport properties in SL. 

In the static regime, the usual method in the investigation of transport properties in SL 
under the influence of external fields is through the observation of static I-V characteristics 
U|. In the realm of ultra-short laser pulse (usually in hundreds of femto-second duration) 
generated coherent regime, what physical quantity can be used to probe the transport prop- 
erties induced by the external fields? This question was partially answered by T. Meier et al 
with the help of the well-known semiconductor Bloch equations || . They found by monitor- 
ing the time-integrated squared THz emission signals that even in the presence of excitonic 
interactions comparable to the miniband widths, the carriers driven by pure intense laser 
fields (ac fields) reveal dynamic localization provided that edFi/hu coincides with the roots 
of the ordinary Bessel function of order zero [0,0. Inspired by their work, we, in this re- 
port, show by numerical simulations that in the dc-ac fields, the detection of time- integrated 
squared THz signals not only can be used to probe the dc-ac induced dynamic localization 
but also be employed to probe the photon-assisted transport in the presence of excitonic 
interactions in the coherent regime. 

We begin with the following semiconductor Bloch equations by including the longitudinal 
external driving fields F(i): 

d e i i i dP(k t) 
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In the above equations, n c (k, t) and n v (k,t) are the populations in the conduction and va- 
lence bands respectively. P(k, t) is the interband polarization. e C] („)(k, t) = e Cj („)(k, t) — 

J2 V(k, k')n C) („)(k / , t) are the electron and hole energies by taking into account of the ex- 
it' 

citonic interaction, while, e c r v )(k,t) are the corresponding energies of the conduction and 
(valence) bands with no Coulomb interaction. fl(k, t) = [J,E(t) + Z)^(k, k')P(k', t) is the 

k' 

renormalized Rabi frequency and V(k, k') is the Coulomb potential in the quasi-momentum 
space. E(t) is the optical field and assumed to be the Gaussian laser pulses ||. F(i) is 
the combined dc-ac fields: F(i) = F + Ficos(o;t). Incoherent dissipative processes are 
phenomenologically described by the last terms in Eqs. (1) and (2). 

If the energy quanta of the optical field E(t) is off resonant with the excitonic resonance, 
the low excitation regime can be produced. In this regime, the perturbation expansion of 
Eqs. (1) and (2) according to the optical field can be performed by the following procedure 
!!: n cM (k,t) = nf} v) (k,t) + nf} v) (k,t) +■ - •; P(k,t) = F«(k,t) + P^(k,t) + .... The 
usual initial condition that SL is in the ground state before the optical pulse is triggered is 
adopted, i.e., n v (k, t = 0) = 1.0, n c (k, t = 0) = 0. Up to the second-order in the optical field, 
we can get the partial differential equations for both the first-order interband polarization 
P^(k,t) and the second-order electron (hole) population density nf) v Jk,t) ||. 

For simplicity, we focus on the 1-D case, which is just the SL growth direction (assumed 
to be z direction). The driving dc-ac field is assumed to be along this direction, and we 
adopt the contact Coulomb potential V8(z — z'), which carries the most important excitonic 
characteristic ||. Under the above assumptions, the partial differential equations for the first 
order P^\k,t) and the second order n^ 2 \k,t) M can be reduced to the following integro- 
differential equations with the help of the accelerated basis in the quasi-momentum k: ( by 
changing k to k — r)(t), and r)(t) is defined as: r)(t) = | J * F{t')dt') [[|. 

j t - % -[e c {k - V (t),t) -e v (k- V (t),t) + V}} pU(k, t) 

= ~^E(t) + V^P^( q ,t)\ - , (3 ) 

9 ^ t] -A^Eit) (4) 



dt 1 h 

where we have introduced the following new notations for convenience: 

P«(M) = PU(k- V (t),t), h% ) (k,t) = n% v) (k-r,{t),t) • (5) 

In Eqs. (3) and (4), we adopt the relaxation approximation by introducing the transverse 
time T\ and longitudinal relaxation time T2 respectively to replace the collision terms. More 



sophisticated description of these dephasing processes calls for Monte Carlo simulation |T0 

As we have mentioned previously, what we want to do is to calculate the time-integrated 
squared THz signal. The THz emission signal in SL can be expressed as S'thz(^) oc d t j( 2 >(t), 
where j^(t) = | X) / dei(k)/dk nf\k,t)dk (i = c,v) is the current due to the non- 

i 

equilibrium distribution of electrons and holes excited by the Gaussian laser pulses. This 
current can be rewritten in the following equivalent form with the help of accelerated basis 
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m=li:j ^ (k - k n{t)) ^\Kt)dk a=c,v). ( 6) 

Solving Eqs. (3) and (4) in the accelerated basis in k space is more convenient than the direct 
integration in the usual quasi-momentum space. In the following numerical simulation, we 
use the similar parameters as those used by T. Meier et al. The combined miniband width 
A = A c + A^ =20 meV (here, we use the tight-binding model as in Ref. [6]); Coulomb 
potential strength 1^=10 meV. The central frequency of the optical field is assumed to be 
located at 2meV below excitonic resonance which satisfies the off-resonant condition [ITT . 
The full width at half maximum of Gaussian laser pulse envelope |i?(t)| 2 is chosen to be 100 
fs. Both transverse and longitudinal relaxation time Ti and T 2 have been set to be 2 ps. 
The energy quanta of the ac field Tiuj is fixed to be 20 meV. To capture the information on 
dynamic localization, we set the ratio of the Stark frequency ujb to the ac field frequency 
uj to be an integer, e.g, n = 1, without loss of generality. We change the ratio eFid/hu 
continuously, and monitor the time-integrated squared THz emission signal through the 
calculation of the following quantity: 

T s <x J dr\S TUz (r)\ 2 . (7) 

The plot of the time-integrated THz signal vs. the ratio eFid/hu has been shown in Fig.l. 
In this figure, it can be clearly seen that when the ratio eFid/hu) scans through 3.95, 6.95, 
10.0, which are approximately the roots of the ordinary Bessel function of the first order J\, 
the time-integrated squared signal falls into valleys. While the peaks in the plot lies in the 
vicinity at those values of eFid/huj that make \Ji{eF 1 d/huj)\ 2 reach the local maxima. From 
the above effect, it is obvious that | Ji(eFid/huj)\ 2 reflects the signal T s profile qualitatively. 
All the above phenomena demonstrate that even in the presence of Coulomb excitonic inter- 
action whose strength is comparable to the miniband widths of SL, dc-ac induced dynamic 
localization still appears. This kind of dynamic localization can be probed through the 
observation of the signal T s just like the case of a pure ac field 0. 

Another aspect of the combined dc-ac fields on the transport of the SL lies in that 
electrons/holes can tunnel through the adjacent wells through the photon emission and 
absorption, when the ratio of the Stark frequency of dc field ujb to the ac field frequency 
uj is an integer ||. This phenomenon is the so-called inverse Bloch oscillators, which was 
found in a recent experimental study by Unterrainer et al in the static regime [[|. The 
coherent regime counterpart of this effect can also be found in our numerical calculation, 
which was shown in Fig. 2. In this figure, we give the plot of the time-integrated squared 
THz emission signal vs. the ratio ujb/uj, and set the ratio eFid/hu to be 2.0. From the 
figure, we can see that the large peaks appear at the location where the ratio ujb/oj is an 
integer. The other distinctive peaks, which appear at fractional ujb/uj can be attributed to 
the dynamical fractional Wannier-Stark ladders theoretically proposed recently P,[T2[ . 



Let us look at the competition effect between the dc-ac induced dynamic localization 
and the photon-assisted transport facilitated by coupling of the combined dc-ac fields. The 
competition between these two effects is shown in Fig.3, where we show the plot of the 
time-integrated squared THz emission signals vs. the ratio ub/uj. In this figure, we set 
the ratio eFid/huj to be around 3.8 which is the first root of the ordinary Bessel function 
of the first order. It can be clearly seen that the signal shows a deep valley when the 
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ratio ujb/uj reaches around the unity. While, the peaks appearing at ujb/uj = 2,3 are still 
present. In other words, the one-photon assisted transport has been suppressed by this 
specific selection of field parameters, and the dynamic localization prevails over the one- 
photon assisted transport. If we shift the ratio from those values which are in the vicinity of 
the roots of the first-order ordinary Bessel function, the one-photon assisted transport can 
revive. This regeneration can be clearly seen from the Fig. 2. We can further suppress the 
two-photon resonance by selecting the ratio ub/u to be 2, and another ratio eF\d/hu) to be 
the roots of the ordinary Bessel function of the order 2, i.e., make J2(eFid/hu) = 0, which 
we don't show for saving space. 

In summary, using the tight-binding model driven by dc-ac fields and perturbatively 
solving the semiconductor Bloch equations in the accelerated quasi-momentum space, we 
found that the time-integrated squared THz emission signal reveals both the dynamic local- 
ization and the multi-photon resonance in the coherent regime. The dynamic localization 
can be fulfilled when the ratio of Stark frequency lob of dc fields to the ac field frequency u 
is an integer n, and the another ratio eF\d/Tiu is located around one of the roots of the n-th 
order ordinary Bessel function J n . The dynamic localization is embodied by suppressing 
the oscillation amplitude of the THz emission signals. Multiphoton assisted transport (res- 
onance) can be identified when time-integrated squared THz emission signals reach peaks 
provided that the ratio lub/uj is an integer, and another ratio eFid/hu is not in the vicinity 
of the roots of the ordinary Bessel function of integral order. Otherwise, the photon-assisted 
transport will be destroyed by the dynamic localization. Experimentally, all the above find- 
ings can be probed by detecting the time-integrated squared THz emission signals from the 
ultra-short laser pulse excited SL driven by combined dc-ac fields. For example, one can use 
the similar expremental arrangement by the UCSB group [|J. 
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FIGURES 



FIG. 1. The plot of time-integrated squared THz emission signal vs. the ratio eFid/hu, 
showing the dynamic localization induced by dc-ac fields. The parameters used in this figure are 
declared in the text. 

FIG. 2. The plot of time-integrated squared THz emission signal vs. the ratio ujb/oj, 
multi-photon resonance can be clearly seen. We choose the ratio eFid/hu) to be 2.0, which is 
well away from the roots of the first-order ordinary Bessel function. 

FIG. 3. The same plot as that of Fig.2, but we choose the ratio eFid/Tiu) to be 3.8, which is 
around the root of the first-order ordinary Bessel functions. 
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